The three-point function of planar quadrangulations 



J. Bouttier and E. Guitter 

Institut de Physique Theorique 

CEA, IPhT, F-91191 Gif-sur-Yvette, France 

CNRS, URA 2306 

j eremie . bouttierOcea . f r 

OO ■ emmanuel.guitter@cea.fr 
O 
O 
(N 

1— i . 

Abstract 

^ I We compute the generating function of random planar quadrangulations with 

three marked vertices at prescribed pairwise distances. In the scaling limit of 



large quadrangulations, this discrete three-point function converges to a simple 
P^. universal scaling function, which is the continuous three-point function of pure 

r-| ! 2D quantum gravity. We give explicit expressions for this universal three- 

^ I point function both in the grand-canonical and canonical ensembles. Various 

limiting regimes are studied when some of the distances become large or small. 
By considering the case where the marked vertices are aligned, we also obtain 
. the probability law for the number of geodesic points, namely vertices that 

[ lie on a geodesic path between two given vertices, and at prescribed distances 

ly-^ ' from these vertices. 

m 
(N 



1. Introduction 



in 
o 

OO . 

Q , 1.1. The problem 



Maps are fundamental objects of discrete mathematics, usually defined as cellular 
^ ' embeddings of graphs into surfaces. They raise many interesting combinatorial and 
■ probabilistic questions. The understanding of their statistical properties is also rele- 
vant to physics, where random maps provide discrete models for fluctuating surfaces, 
for instance in the context of two-dimensional quantum gravity [Q. Many results have 
been obtained over the years for the enumeration of various families of maps, including 
maps carrying spins or particles. Besides Tutte's original approach using some recursive 
decomposition of the maps 0, several new techniques of enumeration were introduced 
over the years, using for instance random matrix integrals |^,^ or, more recently, bi- 
jections with decorated trees like "blossom trees" P-p!0| or "well-labeled trees" [1TT| - |I^ . 
So far, most enumeration results dealt however with global properties of the maps, 
consisting in simply counting these maps or computing global correlations functions for 
some observables, obtained by averaging over the position of the points where these 
observables are measured (see 0). To fully understand the structure of random maps. 
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it is however necessary to have access to local correlations, in which one controls the 



distance between the points where the measure takes place [p!5| , p!6[| . Very little is known 
at this time about local correlations but many results in this direction should be within 
reach by a proper use of the above-mentioned bijection with well-labeled trees. Indeed, 
in this approach, the coding of maps into trees makes explicit reference to the graph 
distance to some origin vertex, thus keeping track in the enumeration of a number of 



graph distances between vertices. A first progress was achieved in Ref. [T^], where the 



so-called canonical two-point function of quadr angulations was computed. This function 
gives the average number of pairs of vertices at prescribed distance from each other in 
the ensemble of quadrangulations (i.e. maps with only tetravalent faces) with a fixed 
number of faces. In particular, for large quadrangulations and at large distances, a 
sensible scaling limit can be obtained where the two-point function converges to some 
universal scaling function, which is the two-point function of pure 2D quantum gravity 
|l^,17|. Another recent progress was made in Ref. [18], where the question of geodesic 



paths between two vertices in random quadrangulations was addressed, and a number 
of results were obtained on the actual dependence of the statistics of geodesies on the 
distance between their endpoints. 

In this paper, we address the question of the three-point function of planar quad- 
rangulations, consisting in enumerating triply-pointed quadrangulations, i.e. quadran- 
gulations with three marked vertices at prescribed pairwise distances. Our main result is 
an explicit formula for the generating function of such triply-pointed quadrangulations 
counted with a weight g per face. Our approach relies on an extension by Miermont of 
the bijection with well-labeled trees which allows to treat the case of multiply-pointed 



quadrangulations |19] via a coding by more general well-labeled maps, which can then be 
enumerated. Considering the scaling limit of large quadrangulations and large distances, 
we then give explicit expressions for the universal scaling form of this three-point func- 
tion, in both grand-canonical and canonical ensembles. In the language of 2D quantum 
gravity, this constitutes the continuous three-point function for the universality class of 
the so-called pure gravity, with the planar topology. We discuss here its main properties. 

The paper is organized as follows: in section 1.2, we start by briefly recalling some 
known results about the two-point function of quadrangulations, and its continuum 
limit. We then present in section 2 our results for the three-point function. We give 
in section 2.1 the explicit formula for the generating function G{di2, (^23, c^si; g) with a 
weight g per face of planar quadrangulations with three distinct vertices, distinguished 
as 1, 2 and 3, at prescribed pairwise distances di2, ^23 and dsi- We then derive in 
section 2.2 its universal scaling form both in the grand-canonical ensemble with a flxed 
"cosmological constant" and in the canonical ensemble of maps with a flxed, but large 
number n of faces. The grand-canonical three-point function is obtained by letting 
g tend to its critical value 1/12 and considering distances scaling as (1/12 — g)~^/^. 
As for the canonical three-point function, it is obtained by letting the size n of the 
quadrangulations tend to inflnity, with distances scaling as n^/^. Various limits of the 
canonical three-point function are discussed in section 2.3, corresponding to the cases 
when one or several of the (rescaled) distances become large or small. Sections 3 and 
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4 are devoted to the actual derivation of our main formula for G'((ii2, ^23, c^ai; fi*)- We 
first recall in section 3.1 the Miermont bijection between, on the one hand, multiply- 
pointed quadr angulations with sources and delays and, on the other hand, well-labeled 
maps. We show in section 3.2 how to use this bijection in the case of triply-pointed 
maps and how to keep track of all pairwise distances by a proper choice of delays. This 
allows us to reduce our enumeration problem to that of counting particular well-labeled 
maps with three faces, with a number of constraints on their labels. The generating 
function for such maps is finally obtained in section 4, where we make use of three main 
building blocks: the already known generating function for well-labeled trees, recalled in 
section 4.1, a new generating function for properly weighted Motzkin paths describing 
chains of such trees corresponding to well-labeled maps with two faces, explicited in 
section 4.2 and another new generating function for so-called Y-diagrams appearing in 
a decomposition of well-labeled maps with three faces, explicited in section 4.3, where 
we finally establish our formula for G((ii2, (^23, d^i', g)- An alternative derivation of the 
above generating function for weighted Motzkin path is presented in Appendix A. We 
discuss in section 4.4 a number of applications of our formula corresponding to the 
so-called (non-universal) "local limit" of large quadrangulations, where we keep the 
distances finite but let the number n of faces tend to infinity. We discuss in details the 
case where one of the three marked vertices is a geodesic point, i.e. lies on a geodesic 
path between the two others. We give in particular the probability law for the number of 
geodesic points at fixed distances from two marked vertices. We finally discuss possible 
extensions of our results and conclude in section 5. 



1.2. The two-point function of quadrangulations 

Before we present our results for the three-point function, it is useful to recall briefly 
some known facts about the two-point function of quadrangulations. This function 
is obtained by considering configurations of doubly-pointed planar quadrangulations, 
i.e. planar quadrangulations with two marked distinct and distinguished vertices. As 
customary, every configuration is counted with a weight g per face and some inverse 
symmetry factor, equal to the inverse of the order of its automorphism group. The 
generating function G{i] g) of such weighted doubly-pointed quadrangulations where 



the two marked points are at distance i from each other was found to be [jT^ 

Ri 



G{i-g) 



log 



logi?i 



with Ri given by fT? 



Rj — R- 



[^ + 3], 



for z > 2 



for i = 1 



i > 1 



XI] 



Here and throughout the paper, we use the shorthand notation: 



1 - 
T 



x" 



X 



(1.2) 



:i.3) 
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The quantities R and x in ( p..2|) are power series in g determined by 

x + l + ^ X + 4+ ^ 



(x + 4 + i) 



2 ' 



R 



a; + 1 + 



1 ' 



:i.4) 



so that we have 
R 

X 



1- y/T^n^ 

l-2Ag- VI - I2g + ^/Q^12g'^ + Qg^l - 12g ^ 



;i.5) 



The quantity G{i; g) constitutes the discrete "grand-canonical" two-point function of 
quadrangulations and corresponds to an ensemble of quadrangulations with a varying 
number of faces, governed by g. We may instead decide to fix the number n of faces 
of the quadrangulation. The corresponding partition function can be obtained via a 
contour integral in g around as: 



Gi^■,g)\,. = ^f^G{^■,g). 



;i.6) 



This quantity constitutes the discrete "canonical" two-point function of quadrangula- 



tions. 



p{D) 



4>(D) 





D 



(a) 



(b) 



Fig. 1: The universal canonical two-point function, given by (a) the prob- 
ability density function p(-D), and (b) the corresponding distribution func- 
tion 

A sensible continuum limit is reached by letting g approach its critical value 1/12 
and jointly taking a large distance z, with the following scaling: 



(7 = — (1-Ae) 
^ 12 ^ ; 



;i.7) 
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where e — > and where A is the so-called "cosmological constant". In this limit, we 
have 

(1.8) 

SO that 

G{i- g) ~ €"^g (d- v^A'/^) with g{D- a) = dDJ'iD; a) = . 

V / sinh (aD) 

(1.9) 

This constitutes the continuous two-point function in the grand-canonical ensemble with 



cosmological constant A [p!5| , p!6|| . Again, we may instead consider the canonical ensemble 
where the number of faces n of the quadrangulations is now fixed to a large value. A 
sensible limit is reached by letting n tend to infinity with the distance between the 
marked vertices scaling as n^/'^. The integral ( p..6|) translates at large n via a saddle 
point estimate into an integral over real values of some parameter ^, upon setting 

^ f 1 + — 1 , i = Dn^/^ . (1.10) 



^ 12 V' ' n , 

Using Eq. ( |1.9[ ) with e = 1/n and A = — and performing a proper normalization 
of G{i;g)\g-n. by the number of doubly-pointed quadrangulations with size n, we then 
obtain a normalized probability density 

r) POO 

piD) = ^ / rf^ e-^' g{D; ^/^3i^) . (1.11) 

This constitutes the continuous canonical two-point function: the quantity p(D)dD mea- 
sures the infinitesimal probability, in the ensemble of large doubly-pointed quadrangula- 
tions, that the two marked vertices be at a rescaled distance in the range [D, D+dD]. Its 
value ( |1 . 1 1| ) is expected to be universal for a large class of maps forming the universality 
class of so-called pure gravity. Finally, another quantity of interest is the (cumulative) 
distribution function 

Q POO 

^(D) = ^ / ee-?V(D; V'=3i|72) (1.12) 

giving the probability that the distance is less than D. The functions p{D) and $(-D) 
are represented in Fig. 1. For small D, we have p{D) ~ (3/7)1^^. 



2. The three-point function 

We now come to our main subject, namely the three-point function of quadrangu- 
lations. In this section, we present our main result, which is an explicit formula for the 
generating function of planar quadrangulations with three marked vertices at prescribed 
pairwise distances. We shall concentrate here on the properties of this three-point func- 
tion, leaving its precise derivation to sections 3 and 4 below. An important application 
is the continuum scaling limit in which the three-point function takes a particularly 
simple universal form. 
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2.1. The three-point function of quadrangulations 

We consider planar quadrangulations with three distinct vertices distinguished as 
1, 2 and 3, and we denote by di2, d23 and dsi their respective pairwise graph distances. 
Prescribing these distances, let G{di2, ^23, <^3i; 9) denote the generating function for the 
number of such triply-pointed quadrangulations, counted with a weight g per face. Note 
that there are no symmetry factors for triply-pointed quadrangulations as the latter have 
no symmetries. The function G{di2, d23, dsi; g) is defined for strictly positive integer 
values of di2, d23 and d^i that furthermore satisfy: 

- the triangular inequalities: di2 < ^23 + '^31 and its cyclic permutations; 

- the condition that di2 + ^23 + d^i is even. 

This latter condition stems from the fact that every planar quadr angulation is bipartite, 
hence the length of any cycle on the map is even. These conditions are best expressed 
through the parametrization 



di2 


= s + t , 




d23 


= t + u , 


(2.1 


d3i 


= u + s , 





easily inverted into 

di2 - d23 + dsi 



t = 



2 


1 


dl2 + d23 - 


d3i 


2 


1 


— di2 + d23 


+ d3i 



(2.2) 




With this parametrization, the above conditions are equivalent to requiring that s, t, 
u be non-negative integers and that at most one of them may vanish. An intuitive 
picture of this parametrization is as follows (see Fig 2): drawing a triangle 123 in the 
plane with sides (ii2, (i23 and d^i, the quantities s, t and u can be viewed as the radii 
of three pairwise tangent circles centered at the points 1, 2, 3. Note that having one 
of the parameters vanishing, say u, corresponds in this picture to having 3 lying on 
the segment 12. On the quadrangulation, this simply means that 3 lies precisely on a 
geodesic path between 1 and 2, so that di2 = d23 + dsi. We shall say in this case that 
the three points are aligned and that 3 is a geodesic point between 1 and 2. 
We find the following expression for the three-point function: 

G{di2, d23, d3i;g) = A,AtA„F(s, t, u; g) , (2.3) 

where As denotes the finite difference operator acting on an arbitrary function f{s) as: 

Asf{s)^f{s)-f{s-l) , (2.4) 

and similarly for and A^, while F{s,t,u; g) has the explicit form: 

F{s,t,u;g) = 

[3li[s + ll[t + ll[u + ll[s + t + u + 3lf (2.5) 

As before, we use the notation ( PT^ ) with a parameter x related to g by Eq. (|1.4|), or 
more explicitly by Eq. (|1.5| ). Note that F{s,t,u;g) is manifestly symmetric in s, t, u, 
which is consistent with G being symmetric in the distances, as expected. 



Relations (|2.3|) and (^]^) hold a priori only in the allowed range of distances, namely 
when s, t, u are non-negative integers and no more than one of them may vanish. In 
the case of aligned points, say when u vanishes, relation (|2.3|) reduces to 



G{s + t, t, s; g) = AsAtF{s, t, 0; g) (2.6) 



as F{s,t, — 1; (7) vanishes identically, while relation (|2.5| ) reduces to 



f , 0; 9) = PI- '^^1- ''^^1- ■ (2.7) 

^ ' 11]. [ii + SL [* + 3], |s + «+lL ^ ' 

It is convenient to slightly extend the domain of definition of G and F as follows: when 
exactly two of the parameters, say t and u, vanish, corresponding to 2 and 3 being 
identical, Eqs. (|2.3| ) and ( |2.5| ) give a result G(s, 0, s; g) = 0, which can still be considered 
as the correct result as we imposed that the three marked vertices have to be distinct. 
Finally, when s = t = u = 0, the above formulas yield the somewhat conventional 
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value G{0,0,0; g) = 1. We shall therefore take the convention that G{d, d,0; g) = 
G{d, 0, d; g) = G{0, d, d; g) = dd,o- With this convention, Eq. (|2.3|) may be inverted into: 



t u 



F{s, t,u;g) = J2J2Yl + ^' + ^' + ^'^ 9) , (2-8) 

s'=0 t'=0 u'=0 

SO that F{s,t,u; g) — 1 can be interpreted as the generating function of quadrangulations 
with three distinct marked vertices at distances di2 = s'+t' , d23 = t'+u' and d^i — u'+s' 
such that s' < s, t' <t and u' < u. In particular, taking the limit s,t,u ^ oo amounts 
to relaxing these constraints, i.e. considering three arbitrary distinct vertices. We have: 

lim F{s,t,u-g) = \+ + f = Ul + {1 - 12g)-'/' 

n>l ^ ^ 

where we identify the coefficient of g"^ as the number of triply-pointed quadrangulations 
with n faces and three distinct marked vertices. 

2.2. Universal continuum limit 

The above formula for the three-point function is easily translated in the continuum 
limit, i.e. when g approaches its critical value 1/12 and all distances are large, with the 
following scaling: 

1 . 

= — 1 - Ae , 

^ 12 ^ ^ ' (2.10) 

di2 = 1^126-^/" , d2s = D23e-'/^ , dsi = i^3ie"'/' , 



with e — > 0. Similarly, the parameters s, t, u defined in Eq. ( |2.1|) have the scaling: 

Se-^/\ t = Te-^'\ s = Ue-^'^, (2.11) 



s 



where 5', T, U are related to D12, -D23 and D31 via relations similar to Eqs. (|2.1|) and 
( ^72|) , namely 

D12 - D2S + Dsi 



D12 


= S + T , 


s 


D23 


= T + U , 


T 


D31 


= U + S , 


U 



2 

1^12 + 1^23 - 1^31 
2 

-1^12 + 1^23 + ^31 



In this limit, we obtain 

G{di2, d2s, dsi;g) ~ e^/^ 2g{Di2, D23, D31; v^A^/^) 



(2.12) 



F{s, t, u; g) ~ e-^/''j^{S, T, U; v^A^ 



/4^ 



(2.13) 
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with finite scafing functions Q and depending on the rescaled distances and on the 
"cosmological constant" A, encoded for later convenience in the parameter 

a = y372AV4 . (2.14) 

Note the factor of 2 in the first fine of Eq. ( p.l3|) which we introduced to compensate the 
fact that this relation holds only for discrete distances satisfying the parity condition 
that di2 + d23 + (iai be even, while G should be taken as zero otherwise. The factor 2 
ensures a posteriori that Q = {2Q + 0)/2 is a correct measure of G on average. Relation 
( P75| ) becomes 

g{D,2, 1^23, Dsi; a) = ^dsdrdu^S, T, U; a) , (2.15) 
while Eq.( |2.5| ) translates immediately into 

3 /' smh{a{S + T + U))smh{aS)smh{aT)smh{aU) Y ifil 
^ ' ' '"^ ~ sinh(a(5 + T))sinh(a(T + [/))sinh(a(t/ + 5)) ; ' ^ ^ ^ 



The above formulas (|2.15|) and ( p.l6|) give the continuous three-point function in the 



grand-canonical ensemble with a varying size of the quadrangulations and a fixed cos- 
mological constant A. This constitutes the grand-canonical three-point function of pure 
2D quantum gravity. 

If we wish instead to work in the canonical ensemble where the number of faces 
n of the quadrangulations is fixed to a large value, and the distances are scaled as 
n^/'^, we can use the above formulas with e = 1/n, and a varying cosmological constant 
A = — As seen for the two-point function, the coefficient of g"" in G or F is obtained 
by a contour integral in g which, asymptotically at large n, yields, via a saddle point 
estimate, an integral over real values of ^, namely: 



G{di2,d23,d3i;g)\g^ ^ / d^^e-^ n'^/^ 20(012, D23, D^i; V-3i^/2) 

-ion /•oo ^ 

F{s, t,u-g)\gr.^— dii e-^ n^/^J^{S, T, U; ^=31^72) , 
mn J_^ 



(2.17) 



where we take the determination V— ir = e ^^snir)nr/4^yj^ ^^j. ^ j-qqI Normalizing these 
asymptotic values, we find the canonical three-point function 

2 

p(Di2,I^23,^3i) = ^ / d^Ce-^'g{Di2,D23,D3i;^/^3ii/^) , (2.18) 

where p(-Di2, -D23, -D31) dD 12 dD23 dDsi is the infinitesimal probability that the pairwise 
rescaled distances between the three marked points in the ensemble of triply-pointed ran- 
dom quadrangulations of fixed large size be respectively in the ranges [-D12, -D12 + dDi2], 
[D23, D23 + dD23], [-D31, -D31 + dDsi]. Eq. ( |2.18| ) constitutes the canonical three-point 
function of pure 2D quantum gravity for the planar topology. 
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Eq. (|2.18| ) was obtained by dividing the first line of Eq. ( p.l7|) by tfie number of 
triply pointed quadrangulations of size n, as given by (|2.9| ): 



3"" f2n\ 12" 



and multiplying by n^/"^ /2 which is the number of allowed triplets of discrete distances 
((ii2, (^235 <^3i) falling into the above range of rescaled distances (note the factor 1/2 
coming from the parity of di2 + (i23 + <^3i)- Similarly, the second line of Eq. ( |2.17|) 
translates into the integrated three-point function 



J-oo 



(2.20) 



corresponding to the probability that the three marked points be at distances S' + T' , 
T + U' and U' + S' with S' < S,T' <T and U' < U. The functions $ and p are 
related through 

p(I)i2, 1^23, Dsi) = ^dsdrdu^iS, T, U) (2.21) 

or conversely 

pS pT p-U 

^{S,T,U)= / dS' dT' dU' 2p{S' + T',T' + U',U' + S') , (2.22) 
Jo Jo Jo 



where the factor 2 may be alternatively understood as the Jacobian in the change 
of variables from (L>i2, -D235 -D31) to (5', T, [/). Finally, Eq. ( |2.y| ) translates into 
limg 7- [/^oo ^{S, T, U) = 1, so that we find the normalization 



dDi2 dD23 dDsi p(Di2, D23, Dsi) = l (2.23) 

x> 

as it should. Here V denotes the domain of positive real values of -D12, -D23 and D^i 
satisfying the triangular inequalities. 

Let us know discuss in more details the properties of the above continuous grand- 
canonical or canonical three-point functions. As a preliminary exercise, let us see how 
to recover the two-point function from the value of the three-point function upon inte- 
grating over the position of one of the points, say 3, keeping the distance D12 fixed to 
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a constant value. We have 



dD23 dDsi g{Di2, D23, Dsi.a) 

D poo poo 

dS dT dU dsdTduJ^{S,T,U-a)5{S + T -D12) 
Jo Jo 

/•oo /"OO 

/ dS dT lim {dsdTJ^{S,T,U;a)) 6{S + T - D12) 
Jo Jo U^oo 



D12 /.i^u^^,Q\.i^u^^.^^> _ q\\\'^ 

dS 18 







sinh(aS') sinh(a(il>i2 — 5")) 
sinh^(ail>i2) 

9 / w ^ sinh(2aDi2l 

'4Di2 + 2Di2Cosh(2aDi2)-3 ^ 



8sinh (Q;i:'i2) V 

(2.24) 

Here T>i2 is a short-hand notation for the domain of aUowed distances -D23, -D31 at 
a fixed value of D12. Comparing this result with the value ( |1.9| ) for the continuous 
two-point function, we see that 

dD23dDsiQ{D,2.D2s.Dsua) = --^d^g{D,2;a) . (2.25) 

Taking a as in Eq. (|2.14| ) for the grand-canonical ensemble, we have —9/{16a^)da: = —Oa 
and we thus find the expected result that the integral of the grand-canonical three-point 
function over the position of one of the points reproduces precisely the grand-canonical 
two-point function for the remaining two points, up to the action of the trivial operator 
—Oa corresponding precisely to the marking of the third point with no prescription on 
its position. 

Similarly in the canonical ensemble, taking a = y^— 3i^/2, we have — 9/(16q!'^)9q: = 
1/(20^5 so that 

f 2 1 

/ dD23dD3ip{Di2,D23,D3i) = ^ d^ ie-^' —d^g{Di2; ^-31^/2) = p{Di2) 

(2.26) 

upon integrating by part over ^, with p{Di2) as in Eq. ( p..ll| ). As expected, we re- 
cover the canonical two-point function as the marginal probability density obtained by 
integrating the canonical three-point function over the position of one of the points. 

In order to visualize the three-point function and in view of the above result, it is 
convenient to first fix one of the distances, say D12 to a constant value and to plot the 
conditional probability density: 

p{D23,D3^\D,2)^ '^'''%'^''f''^ ■ (2.27) 

P{Dl2) 

The quantity p(D23, D3i\Di2)dD23dD3i measures the infinitesimal conditional prob- 
ability that point 3 be at distances from points 1 and 2 respectively in the ranges 
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Fig. 3: Left: plots of the conditional probability density p{D23, -D31I-D12) 
for D12 = 0.8, 1.5 and 3.0, from top to bottom. Right: the corresponding 
contour plots, where the thick lines indicate the boundary of the allowed 
domain V12 for distances, as determined by the triangular inequalities. 
For the first (small enough) two values of -D12, the density is maximal for 
D23 = -D31 ~ 1.5. For the larger value D12 — 3.0, the density becomes 
elongated and squeezed along the boundary line D23 + -D31 = -Di2- 
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[D23,D23 + dD23] and [D^i^D^i + dDsi], given that the distance between 1 and 2 is 
Di2. This conditional probabihty density may then be plotted for various values of D12, 
as illustrated in Fig. 3, where a number of properties of the three-point function can be 
observed. 

First, the three-point function vanishes at the boundary of the domain 1^12, namely 
when the three points are aligned. It also vanishes when the distances from point 3 to 
points 1 and 2 tends to infinity. The probability density is maximal when D23 and D31 
are both equal to a particular value -Dmax, depending on D12. For small enough D12 
(less than ~ 3 or so) , this preferred value is roughly independent of D12 and comparable 
to the maximum of the two-point function p{D), namely D 

max ~ 1.5. For larger D\2'> 

the triangular inequalities impose that -Dmax > -D12/2 so that the value of -Dmax must 
increase. In this case, the density profile becomes squeezed and elongated along the 
boundary line D23 + -D31 = -D12 of the domain 7)12 with, as we shall see below, a width 
of order D^2^^ large D12. 

2.3. Limiting behaviors 

The conditional probability density p(il>23, -D31I-D12) has interesting limiting be- 
haviors whenever D12 becomes small or large. For D12 —>■ 0, taking 5" = aDi2 and 
T = tD 12, we find that 

g{{a+T)D^2, rDi2 + U, aDi2 + U;a) = ^^d^drduJ^{aDi2, rDi2, U; a) 

= ^p;yi QiU, a) + 0{D,,) , 

(2.28) 

where Q{U;a) is the two-point function, as given by Eq. ( |1.9|) . Taking a + r = 1 and 
using the small D12 behavior p{Di2) ~ {?)/7)D\2^ we deduce that 

pi{l-a)Di2 + U,aDi2 + U\Di2) ~ y< p{U) x2la^{l- af {l + 2a{l - a)) (2.29) 

in terms of the canonical two-point function p{U). Writing U = [D23 + -D31 — Di2)/2 
and a = (1 -|- oj)/2 with u = (-D31 — iI>23)/-Di2, this equation may alternatively be 
written as 



,n n in ^ ^ f D23 + D31 - Di2\ . { D31 - D23 
p{D23,D3i\Di2) ~ 7-— X p X i/j 



" V 2 J ^ \ D12 
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with tP{u) = —{l-uY{3-uj^) 



(2.30) 



involving a new function t/^lto) properly normalized to 1 when co varies from —1 to 1. 
The conditional probability density p(L>23, -D31 1-D12) therefore factorizes in the small 
D12 limit into the product of the canonical two-point function in the "longitudinal 
direction" (corresponding to varying values of -D23 + -D31) and of the above probability 
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Fig. 4: The conditional probability density p(L>235 -D31I-D12) for a small 
value of D12, here D12 = 0.05 (a). The same density (b) and its 
contour plot (c), now expressed as a function of the longitudinal vari- 
able U = (-D23 + -D31 — Di2)/2 and a resettled transverse variable u = 
(-D31 — D23)/Di2. As apparent by taking longitudinal and transverse cut 
views along the thick lines in (c), the conditional probability factorizes 
into the product of the two-point function p{U) (red curve in (d)) and the 
probability density '(/'(w) (green curve in (e)). 



density '0 in the "transverse direction" (corresponding to varying D^i — 023)- This 
property is illustrated in Fig. 4 for a value D12 = 0.05. 

Let us now consider the case D12 00, and more precisely the limit when both S 
and T tend simultaneously to infinity (with D12 = S + T). In this limit, we have 

g{Di2, {D12 -S) + U,S + U]a)r^ 24ae-°(2^^2+^^) sinh^ (aU) (3 cosh{aU) + 19 sinh(at/) ) 

(2.31) 

while the grand-canonical two-point function behaves as 

g{Di2; a) ~ 16a^ ^''^^d,^ . (2.32) 
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Fig. 5: The conditional probability density p(il>23, -D31I-D12) for a large 
value of D12, here D12 = 10. (a). The same density (b) and its contour 
plot (c), now expressed as a function of the transverse variable S = (-D12 + 
D31 — D23)/2 and a rescaled longitudinal variable p = (91^12)^/^ (-D23 + 
Dsi —Di2)/2. As apparent by taking longitudinal and transverse cut views 
along the thick lines in (c), the conditional probability factorizes into the 
product of the probability density (^(i/) (red curve in (d)) and a uniform 
density in the 5" direction (green curve in (e)). 

We deduce that 

p{D^2-S + U,S + U\Di2) 

/^(^^^e-^'24ae-«(2i)i2+5f/) ^\^\y^^aU) (3cosh(aC/) + 19sinh(a[/)) 

(2.33) 

with a — 3i^/2 . Note that this result is independent of 5". Both integrals are 
dominated at large -D12 by their saddle point at ^ = {?iD\^^l'^\l2, leading to 



1 .wnr. M/3 J iar^ M/3 ^23 + ^^31-1^12 



p(D23, ^311^12) ~ 777^ X (91^12)'/' ^ (91^12 
^-^^12 \ 

4 

with >^{y) = - sinh(z//2)2 (lle"^'^ - Se"^'") 
3 
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^2.34) 



involving a new function properly normalized to 1 when u varies from to oo. 
In the limit of large D12, the conditional probability density p{D23, -D31I-D12) therefore 
becomes elongated and squeezed along the line D23 + -D31 = -D12, with a profile given by 
the above probability density (p in the longitudinal direction and a uniform probability 
density in the transverse direction. This property is illustrated in Fig. 5 for a value 

— 1/3 

D12 = 10. Note that the extension in the longitudinal direction is of order 5 
announced. 

To end this section, let us finally discuss the value of the three-point function 
whenever all distances become small or large. For S, T and U small (and of the same 
order) , we find the limiting behavior 



9 {STUjS + T + U)f {S^ + T^ + U^ + ST + TU + US) 



(2.35) 



which is a homogeneous function of degree 8 of its arguments. Applying Eq. ( |2.21| ), we 
deduce that p(-Di2, -D23, -D31) is, at small distances, a homogeneous function of degree 5 
of its arguments. Note that for a manifold of fractal dimension (ip, we expect a degree 
{dp — 1) + {dp ~ 2) = 2dF — 3 as the choice of point 2 at a fixed distance from point 1 
will select a manifold of dimension dp — ^ and the choice of point 3 at fixed distances 
from points 1 and 2 will then select a manifold of dimension dp — 2. The value 5 above 
is therefore compatible with the known value dp = 4: for the fractal dimension of large 
quadr angulations. 

Finally, when -D12, -D23 and D^i are large and of the same order, we have 

g{Di2. 1^23, ^31; «) ~ 66 a e-"(^i2+^23+i?3i) (2.36) 
and we find the asymptotic behavior 

p(Di2, 1^23, 1^31) ~ ^(1^12 + 1^23 + I)3i)e-(^)'''(^-+^-+^-)''" . (2.37) 
V6 



3. Triply-pointed quadrangulations and well-labeled maps 

Let us now come to the derivation of our main result, as given by Eqs. ( |2.3| ) and 
( p.5|) . For this derivation, we shall rely on a new bijection discovered recently by Mier- 
mont between multiply-pointed quadrangulations, i.e. quadrangulations with a number, 
say p, of marked vertices, and well-labeled maps with p faces, as will be defined below. 
This bijection extends a well-known bijection obtained by Schaeffer ||Tl|,12| between 



pointed quadrangulations and well-labeled trees (or g-trees) corresponding to the case 
p = 1. In this section, we first recall the Miermont bijection and show how to use 
it to treat the case of planar triply-pointed quadrangulations with prescribed pairwise 
distances between the three marked vertices. We then exploit these results in section 
4 to obtain explicit formulas for various generating functions, leading eventually to 
expressions (|2.3|) and (p75|). 
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3.1. The Miermont bijection 

The Miermont bijection works as follows. We start with a bipartite quadrangulation 
of genus h with n faces, together with a marked p-tuple of distinct vertices, hereafter 
denoted 1,2, and referred to as the sources of the quadrangulation. We shall 

denote by dij, i,j = 1, . . . the distance in the quadrangulation between the sources 
i and j. For the construction below to work, we must impose that dij > 2 for all i j, 
i.e. no two marked vertices are immediate neighbors on the map. 

With each source i, we furthermore associate an integer r^, hereafter referred to 
as the delay of the source i, and which will act as a penalty when measuring distances 
from this source. Delays are subject to the following constraints: 

In -Tj\ < dij, l<i^j<p, ^ ^ 

(3.1) 

Ti — Tj + dij IS even, I < j < p ■ 

Note that delays satisfying these constraints can always be found, for instance by tak- 
ing Ti = (respectively 1) for all sources belonging to the white (respectively black) 
sublattice of the naturally colored bipartite quadrangulation. 

Given a quadrangulation with sources and delays as above, we may associate with 
each vertex v of the quadrangulation a label i{v) defined as 

i{v) = min {Tj+dj{v)) , (3.2) 

where dj(v) denotes the distance in the quadrangulation from the vertex v to the source 
j. The label of any vertex is therefore the minimal value for all sources of a "delayed 
distance" equal, for each source, to the actual graph distance to this source incremented 
by the delay of the source. The first condition in ( |3.1j ) ensures that the label of the 
source i is Ti (with the minimum in ( p.2|) reached only for j = i). Clearly, the labels of 
two adjacent vertices cannot differ by more than 1. Moreover, the parity of Tj + dj{y) 
necessarily changes between neighbors on a bipartite map and, from the second condition 
in ( |3.1|) , this parity is independent of j. This implies that the parity of l{y) must also 
change between neighbors, leading to the crucial property: 

\l{v)-l{v')\ = l (3.3) 

for any pair of adjacent vertices v and v' . The faces of the quadrangulation are therefore 
of two types: simple faces with a cyclic sequence of labels around the face of the form 
£, £ + 1, £ + 2, £ + 1 and confluent faces with a sequence of the form £, £ + 1, i, £ + 1. 
As in Schaeffer's original construction, we associate with each of the n faces of the 
(now labeled) quadrangulation an edge as follows: for each simple face, we select the 
edge of type £ + 2 ^ £ + 1 encountered clockwise around the face (see Fig.6-(a)); for 
each confluent face, we draw a new edge between the two vertices labeled £ + 1 (see 
Fig.6-(b)). 

It was shown by Miermont that the graph spanned by these n edges is a map that 
connects all the original vertices but the sources, which become isolated in the process 
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(a) (b) 

Fig. 6: The two possible types of faces according to ttie labels of ttieir 
incident vertices: (a) simple faces and (b) confluent faces. With each face, 
we associate an edge (thick line) as shown. 

(see Fig. 7 for an example). This map has moreover the same genus h as the original 
quadrangulation, and has p faces, each of them enclosing exactly one of the original 
sources and being referred to as face 1, 2, . . . accordingly. Keeping the labels on the 
vertices of the map, we see that, by construction, these labels vary by at most one along 
any edge of the map. Maps with integer labels satisfying this latter constraint will be 
called well-labeled maps. Finally, it can be shown that, for all the vertices incident to 
face i in the well-labeled map, the minimum in (|3.2|) is reached for j = i (and possibly 



for some other value of j). In particular, all these vertices have a label larger than or 
equal to + 1, and this value is reached by at least one of them, namely: 

min i{v) = r, + 1 . (3.4) 

V incident to face i 



It was shown that the construction above provides a bijection between, on the one 
hand, quadr angulations of genus h with n faces, p distinguished sources and prescribed 
delays {ti} satisfying (|3.1|), and on the other hand, well-labeled maps of genus h with 
n edges and p distinguished faces whose labels satisfy ( |3.4|) . 

To recover the original quadrangulation from the well-labeled map, one can proceed 
as follows: for each face i of the map, we visit all its incident corners successively by 
going counterclockwise around the face. Calling Ti + 1 the minimal label of vertices 
incident to face i, we then connect by an arch each corner with label £ > -f- 1 to the 
first corner labeled i — 1 encountered counterclockwise around the face (see Fig. 8 for 
an example). Corners labeled Ti + 1 are finally connected by arches to an extra vertex 
labeled Ti at the center of the face i. The set of all arches forms the edges of the original 
quadrangulation and the added vertices correspond to its sources. 

An important property of the above coding of multiply-pointed quadrangulations 
by well-labeled maps is that it keeps track of a number of distances in the original map. 
Indeed, as already mentioned, all the vertices incident to face i in the well-labeled map 
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Fig. 7: An example (a) of triply-pointed quadr angulations where the 
marked vertices are indicated by thick circles. Taking delays ti — 0, T2 = 1 
and Ta = 2 as indicated, each vertex v receives its label £{v) as defined by 
(p.2|). Applying (b) the rules of Fig. 6 on all the labeled faces results in a 
well-labeled map (c) with three distinguished faces z = 1, 2, 3. As apparent 
in (c), the minimal label among vertices incident to face i is -|- 1. 

correspond to vertices of the original quadr angulation for which the minimum in ( |3.2|) 
is reached for j = i. This ensures the crucial property: 

di{v) = £{v) — Ti for V incident to face i . (3.5) 

Finally, the above bijection holds for prescribed values of the delays but one may 
as well consider all pairs made of a quadrangulation with p marked vertices and of a set 
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(a) 



(b) 



Fig. 8: The inverse construction leading back from a well-labeled map, 
here with p = 3 faces, to a multiply- (here triply-) pointed quadrangulation. 
As explained in the text, calling Ti + 1 the minimal label of vertices incident 
to face i, an extra vertex (thick circle) with label is added at the center 
of each face i. Each corner with label i is connected (a) by an arch (dashed 
line) to its successor with label i — 1. The set of these arches reconstructs 
the edges of the quadrangulation (b), while the added vertices form the 
sources of the quadrangulation. 

of delays satisfying ( |3.1| ) . In that case, in order to keep a finite number of configurations 
of delays, we must impose an additional condition, say for instance mini=i^...^p — 0. 
The resulting configurations are clearly in bijection with well-labeled maps with p faces 
and arbitrary integer labels such that the minimal label is 1 and the delay associated 
with each source i is recovered via = min„ incident to face i ^{v) — 1. 

3.2. Application to triply-pointed quadrangulations 

Let us now apply the above results to the case of planar quadrangulations with 
three distinct marked vertices, distinguished as 1, 2 and 3. Note first that a planar 
quadrangulation is automatically bipartite, as required by the above construction. We 
shall distinguish two cases of marked vertices: 

(i) the generic case where the three marked vertices are not aligned; 

(ii) the particular case where these three vertices are aligned. 

Let us first concentrate on case (i), leaving the discussion of case (ii) to the end of this 
section. Note first that, since the three points are not aligned, their pairwise distances 
dij, 1 < i j < 3, are all larger than or equal to 2. Indeed, if at least one of these 
distances, say di2, was equal to 1, i.e. if 1 and 2 were immediate neighbors, then, from 
the triangular inequalities, we would have \d3i —d23\ < 1 and, from the bipartite nature 
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(a) (b) (c) (d) 




(e) (f) (g) 

Fig. 10: The seven possible backbones of planar maps with three faces 
distinguished as 1, 2 and 3. It will prove useful to view backbones of type 
(b), (c) and (d) as degenerate limits of the backbone (a) when one of its 
three edges is shrunk to a point. 

of the map, this would imply dsi = + 1 or = (isi + 1. In the first case, 2 would be 
on a geodesic path between 1 and 3 while in the second case, 1 would be on a geodesic 
path between 2 and 3. In both cases, the three vertices would be aligned. 

Having dij > 2 for all 1 < z 7^ j < 3, we may therefore apply Miermont's con- 
struction with p = 3. The resulting well-labeled map is now a planar map with three 
faces that are distinguished. Such maps can be classified as follows: it is convenient 
to view the map as made of a skeleton map decorated by a number of attached tree 
components (see Fig. 9). The skeleton map is obtained by iteratively erasing all the 
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edges of the map which are incident to a univalent vertex until no univalent vertex is 
left. The suppressed edges form tree components which, once attached to the skeleton, 
reproduce the original map. The skeletons themselves can be classified according to 
their backbone map, which is the map obtained by erasing all the bivalent vertices of 
the skeleton and merging their incident edges into a single edge (see Fig. 9). Clearly, 
by construction, the backbones are planar maps with three (distinguished) faces and 
no univalent nor bivalent vertex. There are only a finite number of such maps, all dis- 
played in Fig. 10. This allows to classify all well-labeled maps obtained in the Miermont 
bijection according to their underlying backbone. 

The above classification holds for arbitrary planar maps with three (distinguished) 
faces and makes no reference to labels or delays. For arbitrary fixed values of the delays, 
the maps in the Miermont bijection must be equipped with labels satisfying (|3.4| ) and 
one recovers all quadrangulations with three marked vertices at pairwise distances com- 
patible with these delays by considering all possible well-labeled maps with such labels. 
This requires in particular to consider all possible backbones in the above classification. 

The key point of our derivation that will allow us to keep track of all pairwise 
distances is to impose an extra condition relating the delays at the sources to the pairwise 
distances between these sources. In other words, we choose particular values for the 
delays related to the pairwise distances of three marked vertices on the quadrangulation. 
This particular choice of delays will impose extra conditions on the associated well- 
labeled maps that will restrict the possible choice of backbone and imply additional 
constraints for the labels. More precisely, let us choose for the delays ri, T2 and rs the 
values — s, — t, —u with u as in Eq. ( p.2| ), namely 



—di2 + d23 — dsi 

Ti = — S = , 

2 

r2 = -t= -^12 - ^23 + ds^ ^ ^3_g^ 

di2 — d23 — dsi 
r3 = -u = . 



Note that s, t and u are strictly positive integers since the three points are not aligned, 
so that the delays are strictly negative. We have \ti — T2\ = \d23 — dsi\ < di2 from the 
triangular inequality and again, since we assumed that the points are not aligned, we 
cannot have equality. This property and the same properties under cyclic permutations 
of the indices ensure that the first condition in ( ^.1| ) is satisfied. Moreover, Ti — T2+di2 = 
2t is even and so are the two similar quantities obtained by cyclic permutations of the 
indices, so that the second condition in (|3.1|) is also satisfied. 

Let us now analyze in more details the constraints on the labels of the well-labeled 
maps that result from this particular choice of delays. First, we must impose the general 
condition ( |3.4| ), namely that all labels of vertices incident to face 1 (respectively face 
2 and 3) be larger than or equal to 1 — s (respectively 1 — t and 1 — w), this value 
being reached by at least one of them. In particular, for any vertex v incident to face 
1 (respectively 2 and 3), the quantity i{v) — ri = £{v) + s (respectively £{v) + t and 
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i{v) + u) measures the distance from this vertex to the somce 1 (respectively 2 and 
3). Let us now consider the boundary between, say faces 1 and 2, defined as the set of 
vertices and edges incident to both faces 1 and 2. Note that such boundary is part of 
the skeleton, which is the union of all boundaries for all pairs of faces. For any vertex 
on the boundary between faces 1 and 2, the quantity {i{v) + s) + {i{v) +t) = 2£{v) +di2 
measures the length of a particular path joining the source 1 to the source 2 and passing 
through V. This length must be larger than or equal to the distance di2 between 1 and 
2, which implies that i{v) > 0. The same is true for vertices on the boundary between 
faces 1 and 3 or between 2 and 3. We therefore have the crucial property that the labels 
of vertices on the skeleton are non-negative. 

Now consider a geodesic path on the quadrangulation between sources 1 and 2. 
This path must necessarily intersect the skeleton of the well-labeled map at a vertex vi 
incident to face 1 and at a vertex V2 incident to face 2. We have then di2 = {^{vi) + 
s) + d{vi,V2) + (^(^2) + where d{vi,V2) is the distance between vi and V2 in the 
quadrangulation. As s + t = di2, this implies that i{vi) + i{v2) + d{vi, V2) = 0, and, 
since vi and V2 lie on the skeleton, all these terms are non-negative. We deduce that 
1, V2) = 0, i.e. vi and V2 are identical, and moreover i{vi) = i{v2) = 0. The fact that 
vi = V2 requires that the boundary between faces 1 and 2 be non-empty, which rules 
out well-labeled maps whose backbone is of type (g) in Fig. 10. Similarly, backbones 
of type (e) and (f) are ruled out by considering a geodesic path between 2 and 3 or 
between 1 and 3 respectively. We are therefore left with backbones of the type (a), (b), 

(c) or (d) only. 

In the case of a backbone of type (a), we moreover deduce that there must be 
a label on the boundary between 1 and 2. The same property holds of course for 
the other boundaries too. We end up with the following constraint: each of the three 
boundaries between faces 1 and 2, 2 and 3, and 1 and 3 must contain a vertex v with 
label i{v) = 0. The same constraint holds in the case of a backbone of type (b), (c) or 

(d) . In this case, one of the boundaries is made of a single vertex (which also belongs 
to the two other boundaries and corresponds to the four-valent vertex of the backbone 
or of the skeleton). We deduce that this vertex must have a label 0. 

To summarize the above analysis, the well-labeled maps obtained for our particular 
choice ( |3.6| ) of delays satisfy (see Fig. 11 for an illustration): 



mm 

V incident 
to face 1 



i{v) 



mm 

V incident 
to faces 1 and 2 



mm l[v) 

V incident 
to face 2 



1-t 



mm 

V incident 
to faces 2 and 3 



min i{v) 

V incident 
to face 3 



mm 

V incident 
to faces 3 and 1 



u 



(3.7) 



By the convention that the minimum of the empty set is -|-oo, the last three conditions 
imply that any two faces have a non-empty boundary, hence the backbone is necessarily 
of type (a), (b), (c) or (d). Note also that the first three conditions are simply the 
general requirement ( |3.4| ) in the Miermont bijection, while the last three are specific to 
our choice of delays. For all these conditions to be compatible, it is crucial that s, t and 
u be strictly positive. 
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Fig. 11: A schematic picture of a well-Iabeled map whose backbone is 
of type (a) in the classification of Fig. 10, and the set of constraints on 
its labels for our special choice of delays (|3.6|). This constraints are "face 
constraints" (in squared boxes) for the vertices incident to each face, and 
"boundary constraints" (in ellipses) for the vertices on the boundary be- 
tween two faces (thick lines). The case of maps whose backbone is of type 
(b) , (c) or (d) is obtained from this picture by simply shrinking one of the 
three boundary lines into a single vertex. 



Let us start conversely from a well-labeled map with three (distinguished) faces 
whose labels satisfy ( p.7| ) for some strictly positive integer values of s, t and u. The 
labels satisfy in particular the general condition ( p.4|) for delays ri = — s, T2 = —t and 
r3 = — zi, and applying the inverse construction discussed in section 3.1 will produce 
a triply-pointed quadrangulation with three distinct sources (necessarily at a distance 
larger than or equal to 2) and with associated delays — s, —t and —u. Again a geodesic 
path from the source 1 to the source 2 will necessarily intersect the skeleton of the 
well-labeled map at a vertex vi incident to face 1 and at a vertex V2 incident to face 2, 
so that di2 > (-iyi) + s + ^(^2) + t > s + t. Now, from the fourth condition in (|3.7|) , 
there exists a vertex v at the boundary between faces 1 and 2 with label l{v) = and 
the union of a geodesic path from the source 1 to f and a geodesic path from v to the 
source 2 forms a path from 1 to 2 of length i{v) -|- s -|- i{v) + t = s + t. We deduce that 
the distance between the sources 1 and 2 in the quadrangulation is given by di2 = s + t, 
and similarly ^23 = t + u and d^i = u + s. In particular, as s, t, and u are strictly 
positive, the sources cannot be aligned. 

To conclude, we now have a bijection between, on the one hand, triply-pointed quad- 
rangulations whose three marked vertices are not aligned and have prescribed pairwise 
distances di2, d23 and dsi, and on the other hand, well-labeled maps with three faces 
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i{v) = 



Fig. 12: A schematic picture of a well-labeled map with two faces (a) and 
a marked vertex 3 (big circle) as obtained from the Miermont bijection 
in the case of aligned marked vertices. For our special choice of delays, 
the constraints on labels are two "face constraints" (in squared boxes) for 
the vertices incident to each of the two faces, one "boundary constraint" 
(ellipse) for the vertices on the boundary between the two faces (thick 
line), and a "vertex constraint" that £ = for the vertex 3. Note that this 
picture can be considered as a degenerate limit of Fig. 11 when the two 
boundaries with face 3 reduce to a single vertex. 

whose labels satisfy (\3.1\) for s, t and u related to di2, d23 and dsi via relation (\2.2^ ). 

It remains to discuss the case (ii) of aligned marked vertices. Let us consider the 
case where, say 3 lies on a geodesic path between 1 and 2, i.e. di2 = d23 + (iai, or 
equivalently u = 0. In this CcLSG, cLS the three points are distinct, we have di2 > 2 and 
we may apply the Miermont bijection to our map, now considered as a doubly-pointed 
quadrangulation with sources 1 and 2 only, complemented by an extra marked vertex 
3 lying on a geodesic between 1 and 2. For the two sources, we again choose specific 
delays, namely ri = — s = —dsi and r2 = — t = —d23- We have jfisi — ^23! < (^23+<^3i and 
dsi — (^23 + (^23 + (^31 = 2(i3i is even, so that the delays satisfy (p.l|). The associated well- 
labeled map now has two (distinguished) faces and its skeleton is a simple loop, made of 
all vertices and edges on the boundary of faces 1 and 2 (note that the backbone in this 
case is not strictly speaking a map as it has no vertex). The same argument as before 
shows that, for our particular choice of delays, all labels of vertices on the skeleton are 
non-negative and at least one of them vanishes. Now the vertex 3 is necessarily incident 
to one of the two faces, say face 2. Consider then a geodesic path from 1 to 3. This 
path must meet the boundary between faces 1 and 2 at a vertex v and we have then 
dsi = i{v) -|- s -|- d3{v) where ds{v) is the distance from v to the vertex 3. As s = dsi, 
we deduce that i{v) + ds{v) — and, as both quantity are nonnegative, £{v) = and 
'^3(^) = 0. The vertex 3 necessarily lies on the skeleton and has the minimum label 0. 
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We therefore end up with a weU-labeled map with two (distinguished) faces, with labels 
satisfying (see Fig. 12 for an illustration) 



V incident 
to face 1 



min i{v) = 1 — s , 



V incident 
to face 2 



min i{v) = 1 — t , 



e{v) = , 



(3.8) 



mm 



V incident 
to faces 1 and 2 



and with a marked vertex 3 among those vertices of the skeleton having label 0. 

Starting conversely from a well-labeled map with two faces, and labels satisfying 
( p.8|) for some strictly positive integer values of s and t, and with a marked vertex 3 
on the skeleton with label 0, we can use the inverse construction to recover a quadr an- 
gulation with two sources and a marked vertex. All paths between the sources 1 and 
2 necessarily cross a vertex v on the boundary between faces 1 and 2 so their length 
is larger than {i{v) + s) + {i{v) + t), hence larger than s + t. As the marked vertex 3 
belongs to the boundary and has label 0, it is at distance dsi = s and = t from the 
sources 1 and 2 and the distance between the two sources is exactly di2 = s + t. In 
particular, the three points are aligned with 3 lying between 1 and 2. 

To conclude, we now have a bijection between, on the one hand, triply-pointed 
quadrangulations whose three marked vertices are aligned with, say 3 lying between 1 
and 2, and have prescribed pairwise distances di2, (^23 o-nd dsi = di2 — d23, and on the 
other hand, well-labeled maps with two faces whose labels satisfy ( p.^ for s = d^i and 
t = d23, and with a marked vertex with label on the skeleton. 

Note finally that the case of aligned vertices is not fundamentally different from 
the generic case of non aligned vertices. Indeed, the well-labeled map with two faces 
and a marked vertex of Fig. 12 may be viewed as a degenerate form of the generic case 
of Fig. 11 for which two of the boundaries, say that between faces 1 and 3 and that 
between faces 2 and 3 reduce to a single vertex, so that face 3 in practice disappears 
and its boundary reduces to the single vertex 3. 



4. Enumeration of triply-pointed quadrangulations 

Thanks to the above bijections, enumerating triply-pointed planar quadrangula- 
tions with marked points at prescribed pairwise distances amounts to enumerating well- 
labeled maps having three or two faces, and satisfying ( |3.7| ) or ( |37^ ) respectively. In 
practice, the generic case corresponds to a well-labeled map whose backbone is of type 
(a) in Fig. 10. One then recovers all the other possible cases (backbone of type (b), (c), 
(d) or aligned vertices) by simply allowing one or two of the boundaries to reduce to a 
single vertex. 

In all this section, we shall consider generating functions with a weight g per edge 
of the well-labeled map, corresponding to a weight g per face of the quadrangulation. 
Our approach is as follows. In section 4.1 we recall known results about the enumera- 
tion of well-labeled trees, that is well-labeled planar maps with one face. In section 4.2 
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we consider well-labeled maps with two faces, which correspond to triply-pointed quad- 
rangulations whose marked points are aligned. Considering this simpler situation first 
allows us to introduce some definitions and derive a formula that will be instrumental 
in the next section 4.3, where we consider the general case and derive the expressions 
( ^T3|) and ( p.5|) for the three-point function. Finally section 4.4 is devoted to other ap- 
plications of our formulas, here in the "local limit" for infinitely large quadrangulations 
(in contrast with the continuum limit studied previously). 

4.1. Generating functions for well-labeled trees 

The first building block in our enumeration task is the generating function Ri{g) 
of well-labeled trees planted at a corner with label i and such that all labels are strictly 
positive. With this definition we clearly have Ri = for i < while, for i > 0, Ri 
satisfies the recursion relation 

^ — 9 (Ri-l + Ri + Ri+l) ' ^ ^ 

where by convention, we set Ri\„o = 1 for i > 0. This recursion relation was solved in 



171, with the result 



where the quantities R and x are related to g by Eq. ( |1.4|) , and more precisely Eq. (|1.5|) . 
As already mentioned in the introduction, the quantities log{Ri/Ri-i) for z > 1 and 
log(i?i) for z = 1 can be interpreted as the generating functions for doubly-pointed 
quadrangulations with two marked vertices at distance i from each other (and counted 
with their inverse symmetry factor). This property is a direct consequence of Schaeffer's 
bijection between pointed planar quadrangulations and well-labeled trees. 

Here, from conditions (|3.7|) or (|3.8| ), we will be led to consider well-labeled trees 
planted at a corner with some label i > and such that all labels be larger than, say 
1 — s for some strictly positive integer value of s. The corresponding generating function 
is nothing but Rg+s, as obtained by shifting all labels on the tree by s. Note that, in 
trees counted by Ri+s, the label 1 — s may or may not be reached. 

4-2. Generating functions for well-labeled maps with two faces 

The second building block in our enumeration task corresponds to quadrangulations 
with three marked aligned vertices. By the bijection of section 3, we must actually 
enumerate well-labeled maps with two faces 1 and 2, with labels satisfying (|3.8| ), and 
with a marked vertex on the skeleton having label 0. 

Consider such a well-labeled map. Its skeleton consists of a single loop of arbitrary 
length m, and the sequence of labels along it, read from the marked vertex, can be 
viewed as the sequence of heights in a Motzkin path Ai of length m, i.e. a path going 
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Fig. 13: Illustration of the expression for the generating function 
of well-labeled maps with two faces, corresponding to triply pointed- 
quadr angulations whose marked vertices are aligned. Opening the skele- 
ton (thick edges) at the marked vertex 3, we obtain a chain with non- 
negative labels (£o5 ^1, • • • 5 ^m) and with attached well-labeled trees. Since 
£o = = 0, the sequence of labels along the chain forms a Motzkin path. 
In the enumeration, each vertex with height i in the Motzkin path receives 
a weight Re+s, which counts the possible configurations of the incident 
well-labeled tree in face 1, a weight Ri+t for the incident tree in face 2 and 
a weight g for the subsequent edge on the chain. 

from to in m steps that never dips below height and is made only of +1, or — 1 
steps (see Fig. 13 for an illustration). Let us denote (0 = io^ii, . . . jirn-i^^m = 0) this 
sequence. The whole well-labeled map is uniquely decomposed into the skeleton and its 
attached tree components: we see that there are two well-labeled trees attached at each 
vertex of the skeleton, one surrounded by face 1 with all labels larger than 1 — s, the 
other surrounded by face 2 with all labels larger than 1 — t. The generating function for 
such trees are Rt+s and Ri+t respectively, where i is the label of the attachment vertex 
on the skeleton. This leads to consider the generating function for Motzkin paths with 
such attached trees: 

m — l 

= E E n 9Ri.+s Ri,+t , (4.3) 

Ul^O Motzkin paths of length m j^ — Q 

- A4 = (0=£o, '^m=0) 

where a weight g is attached to each step of the path (corresponding to an edge of the 
skeleton) , in addition to the weight g per edge in the well-labeled trees already counted 
in Ri. Note that there is no weight attached to the last endpoint {k = m) of the path as 
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it should be identified with the first one (/c = 0), so that we do not overcount the trees 
attached there. By convention the path of length m = is counted with a weight 1. 

Conversely, a Motzkin path with such attached trees can be transformed back into 
a well-labeled map with two faces, but this map does not necessarily satisfy ( p?8| ) in 
general, because we have not imposed that the minimum label in face 1 (respectively 
face 2) be exactly 1 — s (respectively 1 — t). Nevertheless, we readily see that these 
global constraints are realized by considering the generating function: 

AgAtXs^t = Xs^t — Xs-i,t — Xs,t-i + Xs-i,t-i (4-4) 

which corresponds to subtracting the contribution of maps where the labels in face 1 are 
larger than 2 — s = 1 — (s — 1) or the labels in face 2 are larger than 2 — t = l — (t — 1). 
In the end, using the bijection of section 3, we find that AgAtXg^t is the generating 
function for triply-pointed quadrangulations whose marked vertices are aligned, with 
pairwise distances s, t and s -\-t. 



It turns out that, with the particular form (|4.2[ ) for i?^, Xg^t has a very simple 



expression. We indeed have the remarkable combinatorial identity: 

^ _ [3]. [^ + 1L [^ + 1]. [^ + ^ + 3]. .... 

[iL[« + 3L[t + 3L[. + t + iL 

which holds for all nonnegative integer values of s and t. This establishes the restricted 
result 



Let us discuss how ( ^75| ) can be proved. First, from its definition (|4.3|) as a sum of 
weighted Motzkin paths, Xg^t satisfies the recursion relation 

Xs,t = 1 + gRsRt Xs,t (1 + gRs+iRt+i (4.6) 

obtained by decomposing the Motzkin path according to its first step. The Motzkin 
path can be of length and receives by convention the weight 1, or it can start with a 
or +1 step yielding a first factor gRgRt- If it starts with a step, the rest of the path 
is again a Motzkin path with generating function X^^f If is starts with with a step +1, 
the rest of the path is now a path from height 1 to height which, upon considering the 
first return at height 0, can be decomposed into a path from 1 to 1 with heights larger 
than or equal to 1, a step down from 1 to and a Motzkin path from to 0. These 
three components have respective weights Xs+i^t+i, gRs+iRt+i and Xs^t- Gathering all 
the possibilities above leads to Eq. ([4.6|), which uniquely determines Xg^t for all s,t > 
as power series in g, with Xg^t = 1 + 0{g). Now it is a straightforward exercise to 
check that the form ( ^.5| ) of Xg^t does indeed satisfy ( |4.6| ) for Ri given by ( [4.2| ) with the 
relation ( p..4|) , and has the expansion Xg^t = 1 + C'(fi')- This completes the proof. 

The above argument does not allow to derive the explicit form (|4.5|) ab initio from 
its definition ( [4.3| ) but only checks a posteriori that this form matches the definition 
of Xg^f A more constructive and historical proof of the formula ( |4.5| ) is presented in 
Appendix A where we give two methods for enumerating so-called "quadrangulations 
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with a geodesic boundary" of length 2(s + t). Equating the two results precisely yields 
the expression ( ^?5| ) for Xg^t- 

Yet another proof is provided by the discussion in section 4.4 below, where it will 
be shown that AgAtlogXg^t is the generating function for well-labeled maps with two 
faces satisfying ( |3.8| ), but without a marked vertex on the skeleton (this induces non- 
trivial symmetry factors). By the Miermont bijection, with the particular choice of 
delays ri = —s and r2 = — t, these correspond to doubly-pointed quadrangulations 
where the marked vertices are at distance s + 1. But the generating function for these 
doubly-pointed quadrangulations is given by Eq.(|l.lj) with i = s + 1, and we arrive at 
the equality: 

log (^) = A.A, logX., = log ( ^-'■^'-■■'- ) . (4.7) 

It is now an easy exercise to derive ( ^4.5| ) from this equality using ( |4.2| ) and the initial 
values Xs,o = Xo^t = 1- 

Finally, let us remark that: 

lim Xs,t = 1 + X + = . (4.8) 

Thus the basic quantity x can be interpreted as a generating function for Motzkin paths 
with attached well-labeled trees with no lower bound on the labels. 

4-3. Generating functions for well-labeled maps with three faces 

We are now ready to enumerate triply-pointed quadrangulations in the generic 
case, corresponding to well-labeled maps with three faces 1, 2, 3 satisfying (|3.7|). Let us 
consider such a well-labeled map. As discussed in section 3, its backbone is necessarily 
of one of the types (a), (b), (c) or (d) illustrated in Fig. 10. 

We first consider the generic type (a), and denote by v and v' the two trivalent 
vertices of the backbone, say with v such that faces 1, 2,3 appear in counterclockwise 
order around it. Then the skeleton consists of three chains connecting v to v', corre- 
sponding to the three boundaries 1 — 2, 2 — 3 and 3 — 1. By ( p.7|) , the minimal label 
on each boundary must be 0. Let us then denote by Vi2 (resp. v'12) the first vertex 
with label encountered when following the boundary 1 — 2 starting from v (resp. v'). 
Similarly we define V23, t'23, vsi and ^3^. We can now perform a decomposition of the 
well-labeled map by cutting it at the vertices V12, f 12? "^^23, V31 and v'^i, as illustrated 
on Fig. 14. More precisely for each cut vertex, the cutting line is drawn between the two 
corners following immediately the two incident skeleton edges, going counterclockwise 
around the vertex. With this convention, each half of the cut vertex has one incident 
skeleton edge, and one attached well-labeled tree. In general the skeleton is split into 
five components: 

- three (linear) chains connecting V12 to f^2 5 "^23 to ^235 ^si to v'^i respectively, 

- two Y-diagrams, made of three chains connecting a central vertex [v or v') to 
distinct endpoints (f 12, 'i^23, "f^si or '^'^2; "^23; "^31)- 



30 




Fig. 14: The decomposition of a generic well-labeled map with three faces 
into chains and Y-diagrams. We represented the skeleton of the map by 
thick lines and schematized the attached trees by blobs. The skeleton is 
the union of the three boundaries 1 — 2, 2 — 3 and 3 — 1 which meet at 
vertices v and v'. On each boundary, we select the vertices with label 
closest to V and v' (denoted fi2, fi25 "^^23, '^'235 "^31 ^3^^ with obvious 
conventions). We cut the map at these selected vertices as shown by the 
dashed lines. We end up with five pieces in general, namely three chains 
(in red) and two Y-diagrams (in blue). Taking into account the constraints 
on the labels of the attached trees (see text), this decomposition translates 
into the relation F{s, t, u; g) = Xs^tXt^uXu,s(Ys,t,uY ■ 



There are a number of particular cases. If there is a unique vertex with label on one 
of the boundaries, say 1 — 2, then f 12 = f 12? and the corresponding linear chain can be 
seen as an isolated vertex without attached trees, i.e. it is trivial. If, say, v has label 0, 
then V — V12 — V23 = f3i, and the corresponding Y-diagram can also be considered as 
trivial, made of an isolated vertex without attached trees (using our convention to cut 
within the corners following immediately the incident skeleton edges counterclockwise, 
the three trees attached to v go in different chains). For a backbone of type (a), we 
have V ^ v', hence it is not possible that both two Y-diagrams and one of the chains 
be simultaneously trivial. Such situation however occurs when we consider degenerate 
backbones types (b), (c), (d), in which case both two Y-diagrams and exactly one of 
the chains are trivial, or when we consider well-labeled maps with two faces, for which 
both two Y-diagrams and exactly two of the chains are trivial. It is then convenient to 
also include the case of the fully trivial map reduced to a single vertex, when both two 
Y-diagrams and three chains are trivial. To summarize, gathering all possible cases, an 
arbitrary triply-pointed quadrangulation is in correspondence with a well-labeled map 
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that itself can be decomposed into a combination of three chains and two Y-diagrams 
(possibly reduced to isolated vertices). 

Let us now investigate the constraints on the labels of attached well-labeled trees 
inherited from the first line of (|3.7D . It is convenient to slightly relax these conditions 
and consider well-labeled maps satisfying: 

min i{v) > 1 — s , min i{v) > 1 — t , min i{v) > 1 — u . (4-9) 

V incident v incident v incident 

to face 1 to face 2 to face 3 

We shall call F{s,t,u;g) the generating function of such maps. Then the generating 
function for maps satisfying (|3.7| ) is simply given by AsA^A^F. With the relaxed 
conditions ( [4.9|) , the constraints on labels simply factorize into individual constraints 
for each attached well-labeled tree. 

The case of chains is easy: the sequence of labels along the chain between Vi2 and 
v'i2j say, can be seen as a Motzkin path, and all attached trees previously surrounded 
by face 1 (resp. 2) have labels larger than 1 — s (resp. 1 — t). Upon closing the chain, 
this is the same object as that considered in the previous section since, for non-trivial 
chains, there is only one attached tree at the endpoint Vi2 (resp. ^12)5 with labels larger 
than 1 — s (resp. 1 — t). Therefore the generating function for chains with attached 
well-labeled trees satisfying these constraints is nothing but Xg^t, including the trivial 
chain with weight 1. By the same reasoning, upon a common cyclic permutation of 
1,2,3 and s,t,u, we find that the generating function for all possible chains between 
V23 and V22. (^6sp. -^si and v'^-^) is Xt,u (resp. X„,s). 

We next need to consider the new situation of a Y-diagram, say the one containing 
V. As said above, it is made of three chains (branches) connecting the central vertex v, 
to the endpoints t'i2, ^23? ^31- If "i^ has label then the Y-diagram is trivial, otherwise 
every branch has a non-zero length, since i{vi2) = ^(^23) = ^(^31) = 0. On each branch, 
the sequence of labels can be viewed as a path going from height i{v) to 0, made only of 

or — 1 steps, and not reaching height before the last endpoint. On the branch to 
V12, every intermediate vertex has two attached trees, and by ( [1.9| ) the one (within face 
1) has labels larger than 1 — s and the other (within face 2) has labels larger than 1 — t. 
The final vertex V12 has only one attached tree with labels larger than 1 — t. Similar 
properties are found for the other branches, upon a simultaneous cyclic permutation of 
1,2,3 and s,t,u. Finally v, being incident to all three faces, has three attached trees 
with labels respectively larger than 1 — s, 1 — t and 1 — u. Let us denote by Y^^t^u 
the generating function for Y-diagrams satisfying such constraints. This function is 
determined by the recursive equation: 

Ys,t,u = 1 + 9^RsRs+iRtRt+iRuRu+iXs+i,t+iXt+i,u+iXu+i,s+iYs+i,t+i,u+i (4.10) 

which is derived in the same spirit as Eq. (|4.6|) for Xg^t- Indeed, let us consider the 
label of f . If i{v) = then the Y-diagram is reduced to an isolated vertex with weight 
1. Otherwise i{v) > 1, and we cut the Y-diagram at each first vertex with label 1 
encountered when following each branch starting from v (using the same procedure as 
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above for dispatching attached trees). We then obtain four pieces : another Y-diagram 
(which is trivial iff i{v) = 1) and three chains. The Y-diagram has labels 1 on its 
endpoints, otherwise it satisfies the same constraints regarding the labels on its attached 
trees. By decreasing all labels by 1, we find that such Y-diagrams are enumerated by 
Ys+i^t+i,u+i- The chains can be seen as paths from 1 to 0, that do not reach height 
before the endpoint. By cutting out the last step and decreasing all labels by 1, we 
obtain Motzkin paths enumerated by Xs+i^t+i, Xt+i,u+i, -^u+i,s+i, for the respective 
endpoints fi2, f23 and v^i. Finally the last steps contribute with respective weights 



gRs+iRt, gRt+iRu and gRu+iRs- Collecting all contributions, Eq.( |4.10[ ) follows. 
We found the remarkably simple form for Ys^t,u'- 

Y = [^ + 3].[t + 3]Jt. + 3],[. + t + z. + 3], 

[3l[s + t + 3l[t + u + 3l[u + s + 3l ^ 

which can be readily checked by substituting into ( f4.10|) , and noting that it is the unique 
solution satisfying Ys^t,u = 1 + C>{g). However we do not have a more "combinatorial" 
derivation for this formula, similar to those mentioned in the previous section for Xs^t- 
Note that the form ( |4.2| ) for Ri itself still lacks such a combinatorial explanation. Al- 
ternatively, we may also write Ys^t,u as a sum: 

oo 

^=0 (4.12) 

[s + 11 [s + 21 [tl [t + 3l[2i + s + t + 31 



s,t 



where Xi^s,t satisfies Xi^s,t = 9Rs+iRtXs+i,t+iXi-i^s+i,t+i, and corresponds to the 
generating function for branches from v to f 12 where v has a prescribed label i. 

So far we have derived the generating function for Y-diagrams containing v, however 
we easily see that Y-diagrams containing v' have the same generating function Ys^t,u 
(which is symmetric in s,t,u). Combining the expressions for X and Y, we arrive at 
the following expression for F{s,t,u; g): 

\2 



F(s, t, u; g) = Xs,tXt,uXu,s (^s,t,u)' 



mA[s + i]At + MAu + MAs + t + u + 3i 



\2 



[l]'[^ + ^ + l]J^ + ^ + 3]Jt + ^ + l]Jt + ^ + 3L [u + s + ll[u + s + 3l 

(4.13) 

and we readily recognize ( p.5|) . As mentioned above, the desired generating function 
for maps satisfying (|3.7|) is AtA^F(s, i, w; (7). By the discussion of section 3, this is 



precisely the three-point function G{di2, d23, dsi, g). This completes the proof of (|2.3|) . 

4.4- Local limit for large quadrangulations and statistics of geodesic points 

As a final application of our formulas for triply-pointed quadrangulations, we can 
consider the "local limit" of large quadrangulations, obtained by considering the canon- 
ical ensemble of quadrangulations with n faces and letting n tend to 00, keeping the 
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distances between the marked vertices finite. As opposed to the scaling limit discussed 
in section 2.2, results in this local limit are expected to be non-universal and specific 
to quadr angulations, and would be different for other classes of maps. Again, we can 
extract the term (7" of the various generating functions by contour integrals in g of 
the type of Eq. ( p..6|) . In the local limit of large quadrangulations, these integrals can 
be evaluated exactly by saddle point estimates (see for instance Ref. [jT^ for a general 
discussion of this technique). By proper normalizations, our enumeration results then 
translate directly into average properties in canonical ensembles of simply- or doubly- 
pointed large (meaning strictly speaking with fixed size n 00) quadrangulations, i.e. 
quadrangulations with one or two marked vertices, here referred to as sources (see again 
Ref. fl^ for a discussion on these ensembles). 



For instance, from the explicit form ( p.5| ), we can extract the average number 
(c((ii2, (i23, c^si)) of couples of vertices (2,3) at finite distance d23 from each other and 
at respective distances di2 and dsi from the source (denoted 1) in the ensemble of large 
simply-pointed quadrangulations. We find the following formula: 



f{s,t,u] 



(c((ii2, ^23, f^si)) = AsAfAufis, t, u) , where 

_ _9 {{l + s){l + t){l + u){3 + s + t + u)f 

~ 140 (1 + s + t)(3 + s + t){l + t + u){3 + t + u){l + U + s){3 + u + s) 

X (29 + 20(s + t + u) + 5(s^ + t'^ + u'^ + st + tu + us)) 
X {{st + tu + us + stu){4 + s + t + u) — stu) 



(4.14) 

and s, t, u are related to (ii2, ^23, f^si via (|2.2|). 

Taking u — Q, this quantity reduces to AsAt/(s,t, 0) with 

9 (l + 8)(l + t)(3 + . + t) 
' ' ^ 140(3 + s)(3 + t)(l + s + t) (4.15) 

X st (29 + 20(s + t) + 5(s2 + + st)) {A + s + t) 

and measures, in the ensemble of large simply-pointed quadrangulations, the average 
number of pairs made of a first vertex at distance s -\-t from the source and a second 
vertex lying on a geodesic path between the source and the first vertex, at distance 
s form the former. Upon dividing by the known average number A^s+t of vertices at 
distance s-^t from the source in large simply-pointed quadrangulations (see for instance 



Ref. IjT^ for the value of Ns+t), this gives the average number 

(4.16) 



(c(s)),+t = -3— A,AJ(s,t,0) 



with Ns+t = ^{s + t + l)(5(s + t)2 + 10(s + t) + 2) , 
35 

of geodesic points, i.e. vertices lying on a geodesic path between the two sources, and at 
distance s from the first one, in the ensemble of large doubly-pointed quadrangulations 
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Fig. 15: The "profile of geodesic points" {c{s))d, measuring the average 
number c of vertices lying on a geodesic path between the two sources, 
and at distance s from the first source, in the ensemble of large doubly- 
pointed quadr angulations with the two sources at distance d from each 
other. This profile is represented in (a) versus s/d for d = 10, 100 and 
1000 (from bottom to top), for all allowed values of s. The edge of the 
profile corresponding to values s = 1, . . . , 10 is represented versus s in (b) 
for d — 10, 20 and 50 (lower values, ordered from bottom to top) and in 
the limit d oo (top values in red) . 

where the two sources are at distance s + t from each other. Fixing this distance to 
a constant value s + t — d and letting s vary, the "profile of geodesic points" {c{s))d 
versus s/d is represented in Fig. 15-(a) for d = 10, 100 and 1000. When t becomes large, 
keeping s finite, we have 

9 sil + s). 



and thus 



A,AJ(s,t,0) 



28 3 + s 



9 s{5 + s) 



7{3 + s){2 + s) 



(4.17) 



(4.18) 



which compared to the average number {3/7)t^ of vertices at a large distance (s + t) ~ t, 
gives an average number 

, , „ 3s(5 + s) , , 

C « oo = ^ 4.19 

(3 + s){2 + s) 

of geodesic points at finite distance s from the first source when the second source is 
far apart. This limiting form is represented in Fig. 15-(b), where it is compared with 
the profile at finite values d = 10, 20 and 50 of the distance between the two sources. 
In particular, far from the first source, i.e. when s becomes large, Eq. ( |4.19| ) gives an 
average number of 3 geodesic points far away from both sources, as apparent in Fig. 15- 
(a) and (b). 
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Beside its average {c{s))s+t, we have access to the fuU probabihty law for the number 
c of geodesic points. Indeed, consider triply-pointed quadr angulations where the three 
sources are aligned, with the third one lying on a geodesic path between the first two. 
As discussed in section 3.2, these are in bijection with well-labeled maps of the type of 
Fig. 12, with two faces and labels satisfying (|3.8| ). As already discussed, these maps are 
counted by AsAtXg^t, where Xs^t accounts for the closed chain of vertices forming the 
boundary of the two faces. Now, for a fixed well-labeled map, the third marked vertex 
(vertex 3 in Fig. 12) can sit on any of the vertices of the boundary with label and 
the number of geodesic points is given by the number of such vertices. Maps whose 
boundary has exactly c such vertices are simply enumerated by AgAtXl'^^, where 

^i^ = -('^^V. (4.20) 



c V Xs^t / 

Indeed, the quantity {Xs,t-l)/Xs,t = gRsRt{'i-+9Rs+iRt+iXs+i,t+i) is the generating 
function for weighted Motzkin paths having all intermediate heights strictly above and 

(c) 

XgJ simply counts cyclic sequences of exactly c such paths. In particular, we have the 
consistency relation 



J2cX 

c>l 



- X,,t (4.21) 



since in Xg^t, one of the c vertices with label is marked on the chain. We also have 
the normalization 

5^xiJ = logX,,, (4.22) 

C>1 

which allows to recover the generating function for doubly pointed maps with two 
marked vertices at distance s + t via 

A A 1 ^ 1 fXs,tXs-l,t-l\ , f [s + t + S]Js+t]l \ ^ [ Rs+t 

AsAt\ogXs,t = \og[ — ^ =log ^ =log' 



Xs,t-iXs-i,tJ \[s + t-l]^[s + t + 2\l.J \Rs+t-l, 

(4.23) 

As mentioned in section 4.2, this identity can be used to actually derive formula (|4.5| ). 
Combining ( [4.20| ) and ([4.23| ), we deduce that, in the canonical ensemble of doubly- 



pointed quadrangulations with n faces whose two sources are at distance s + t from each 
other, the probability Ps+t;n{cT s) that there be exactly c vertices lying on geodesies 
between the two sources at distance s from the first one reads 

A^AfX^ t 

f-'-<<=-^'=iog(ij,+,/ij:+,_,)i,» ■ 

Again we can consider the local limit of large quadrangulations by sending n ^ oo and 
keeping s and t finite, leading, after extracting the g"^ terms by a saddle point technique. 
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to a probability 



1 ^ ^ /(s,t,0) ' 



with Ag^t = 3 



(4.25) 

(g + l)(t + l)(g + t + 3) 
(s + 3)(t + 3)(s + t + l) 



and with /(s, t, 0) and A^s+t as in ( [4.15| ) and ([4.16| ). In particular, when t becomes large, 
we find a probability 

that there be c geodesic points at distance s from the first source if the two sources 
are far apart. Finally, far away from the first source, i.e. when s becomes large, this 
probability reduces to 

Pooic) = l(l] , (4.27) 



2 V3^ 

which is the probability law for the number c of geodesic points at any fixed distance, 
far away from the two sources. 



5. Discussion 

In this paper, we computed the generating function for triply-pointed planar quad- 
rangulations with three marked vertices at prescribed pairwise distances. We then 
derived its universal scaling form and analyzed its behavior in various limiting regimes. 
The main ingredient in our derivation is the Miermont bijection between triply-pointed 
planar quadrangulations with sources and delays and well-labeled planar maps with 
three faces. To keep track of all distances between the marked points, we had to sup- 
plement this bijection with a particular choice of delays, resulting in a more restricted 
set of well-labeled maps amenable to a direct enumeration. The combinatorial building 
blocks in this enumeration were the already known generating function Ri for well- 
labeled trees and new generating functions Xg^t and Ys^t,u describing configurations of 
chains or Y-diagrams of such trees. These new functions have very simple forms ( [4.5|) 
and (|4.11| ), the first one being also used to address the statistics of geodesic points in 
doubly-pointed quadrangulations. We believe that these formulas will be useful to other 
enumeration problems related to well-labeled maps, and consequently to quadrangula- 
tions, including maps with higher genus. In this respect, is it tempting to view Xg^t as 
a propagator and Ys^t,u as a vertex, whose combination into Feynman diagrams builds 
general well-labeled maps with a number of constraints in their labels. 

A natural question is of course that of the p-point function for p > 3, which would 
require to compute the generating function for multiply-pointed quadrangulations with 
p sources at prescribed pairwise distances. Even if the Miermont bijection works in 
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this case with arbitrary delays, it does not seem possible in general to keep track of all 
pairwise distances by a proper choice of delays as we did for p = 3. Indeed, there are 
only p free values for the delays but p{p — l)/2 pairwise distances, so that, when p > 3, 
we cannot encode all pairwise distances in the delays. More precisely, we expect that 
our construction will work only if the p marked vertices have pairwise distances which 
can be realized as distances between the centers of p pairwise tangent hyper-spheres in 
p — 1 dimensions. In other words, this requires that we can write dij = Si + Sj for some 
set of no n- negative integers Si, i ~ 1, . . . This is clearly not a generic situation but 
it includes for instance the case of regular simplices, for which all pairwise distances are 
the same. Multiply-pointed quadrangulations whose p sources have prescribed pairwise 
distances of the above restricted form can be enumerated explicitly along the same 
lines as in the case p = 3. For instance, taking (81,82,53,84) = in the case 

p = 4, one finds a generating function A^A^ A^A„F(s, t, m, f ; (7) with F{s,t,u,v; g) = 

Finally, another extension concerns classes of maps more general than quadrangu- 
lations. The discrete two-point function is known for instance for bipartite maps with 
prescribed face degrees |]T^ and has a similar structure as that of Eqs. (|1.1| ) and (|1.2| ) , 



but now with [i] taking a more general explicit form. We expect that the Miermont 



bijection naturally extends to these maps, using the mobile rules of Ref. [|T^ . We further 



expect that many of our arguments hold in this case. In particular, the derivation of 



Xg^t in Appendix A should extend to this situation and lead to the same formula (|4.5|) , 
now with the modified [i]^. As for Ys^t,u, we had no combinatorial derivation but we 
may still speculate that the form ( [4.11| ) remains unchanged. 
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RAGE European network, MRTN-CT-2004-5616 and from the Programme d'Action 
Integree J. Verne "Physical applications of random graph theory". 



Appendix A. Derivation of formula (|4.5|) for Xs^t 

We give here a derivation of the explicit formula ( |4.5D for Xg^t based on the enumer- 
ation of quadrangulations with a geodesic boundary. Such objects were introduced and 
enumerated in Ref.]!^] in the context of a general study of geodesic paths in quadran- 



gulations. They can be defined as follows: a quadrangulation with a geodesic boundary 
of length 2i [i > 2) is a planar map with a marked face of degree 2i and with all other 
faces of degree four (squares), and with two marked (and distinguished) vertices incident 
to the marked face at distance i from each other in the map. This last requirement is 
equivalent to demanding that the boundary of the marked face be made of two paths 
of length i joining the marked vertices, which are moreover geodesic paths in the map. 
These paths may possibly meet at common vertices (necessarily at the same distance 
from the marked vertices) or even along common edges. It is convenient, when drawing 
this map in the plane, to choose the marked face as the external face. Ignoring this face 
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Fig. 16: A schematic picture (a) of a quadr angulation with a geodesic 
boundary of length 2i. The underlying quadrangulation is symbolized by 
the grey background. Its boundary is made of two paths of length i joining 
two marked vertices (thick circles) which are geodesic paths in the whole 
quadrangulation. The configuration can be transformed into a true quad- 
rangulation without boundary by adding winding edges as in (b) or (c). 
Such quadrangulation can be considered as a simply-pointed quadrangula- 
tion (b) whose source is the top marked vertex to which we assign a delay 
0. Performing the Miermont bijection creates a well-labeled tree whose 
structure is displayed in (d). The quadrangulation can be considered al- 
ternatively as a doubly-pointed quadrangulation (c) whose sources are the 
two marked vertices to which we assign delays —s and — t, (with s + t = i). 
Performing the Miermont bijection now creates a well-labeled map with 
two faces whose structure is displayed in (e). 
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results in a quadrangulation having a boundary made of the two geodesic paths above 
joining the marked vertices (see Fig. 16- (a) for an iUustration) . 

To enumerate these maps, we may transform them into true planar quadrangula- 
tions by adding edges in the external face linking the two boundaries of the quadran- 
gulation. We shall refer to these boundaries as the left and right boundary by viewing 
the quadrangulation with the first marked vertex at the top. We then connect by an 
edge winding around the quadrangulation each vertex of the right boundary at distance 
k — 2, . . . , z — 1 from the first marked vertex to the vertex of the left boundary at 
distance k — 1 form the first marked vertex (see Fig. 16-(b) and (c)). As a result, the 
external face is divided into i — 1 squares. Note that the added edges do not modify the 
distances to the two marked vertices of any vertex in the map. In particular, the two 
boundaries remain geodesic paths. 

We then can use the Miermont bijection to transform this quadrangulation into a 
well-labeled map. Let us do this in two ways: 
(i) by considering the map as a pointed map with only one source (the first marked 

vertex) and delay ri = (this corresponds to the original Schaeffer construction). 

This construction will create a particular type of well-labeled maps with one face. 



i.e. a well-labeled tree referred to as a "spine tree" in Ref. . 
(ii) by considering the map as a doubly-pointed map with two sources (the two marked 
vertices) and delays ti = —s and T2 = —t for arbitrary strictly positive values of s 
and t satisfying s + t = i. This will produce a particular type of well-labeled map 
with two faces. 

In the first case, the label of the vertices of the two boundaries are their distance to 
the first marked vertex, ranging from to a maximal value i for the second marked 
vertex. Applying the construction of Fig. 6, each of the added squares selects an edge 
of the right boundary (see Fig. 16-(d)). All the edges of this boundary but the first 
one are selected, creating a chain of length i — 1 whose vertices have labels 1, 2, . . . , i. 
The rest of the quadrangulation gives rise to tree components that are attached only to 
the left side of the chain. All these tree components are well-labeled trees with labels 
larger than ri -|- 1 = 1, i.e. with strictly positive labels. Note that the minimal label 1 is 
automatically attained at the first vertex of the chain. The generating function for the 
obtained well-labeled map with the above structure, counted with a weight g per edge, 
is simply: 

Thanks to the Schaeffer bijection, this is also the generating function for our quadran- 
gulations with a geodesic boundary (up to a factor g^~^ for the spurious added squares). 

In the second approach, the labels of the vertices on the boundaries are, when going 
away from the first marked vertex towards the second marked vertex, — s, 1— s, 2— s, . . . , 
for the first s -|- 1 vertices and then —1, —2, . . . ,—t for the last t. Applying again the 
construction of Fig. 6, the first s — 1 added squares select edges of the right boundary 
while the last t — 1 ones select edges of the left boundary (see Fig. 16-(e)). The first s 



40 



edges of right boundary but the first one are selected, creating a chain of length s — 1 
whose vertices have labels 1 — s,2 — s,...,0, while the last t edges of left boundary 
but the last one are selected, creating a chain of length t — 1 whose vertices have labels 
0,— — t. Finally, the s-th added square is a confluent face that gives rise to a 
winding edge connecting the two vertices labeled on the two boundaries. The rest 
of the quadrangulation creates a path joining these two vertices with label 0, and a 
number of tree components attached to both side of this path, to the left side of the 
chain of length s — 1 and to the right side of the chain of length t — 1 (see Fig. 16-(e)). 
All these tree components are well-labeled trees with labels larger than ri + 1 = 1 — s 
or r2 + 1 = 1 — t according to which face they lie in (see Fig. 16). Note again that the 
minimal label 1 — s or 1 — t is automatically attained on the chains. The generating 
function for the obtained labeled maps with the above structure is now 

„s+t-l TT P TT P s+t-l T?s+t + iMx + '^]x y / A oN 

' il " PLl^ + iLPLit-n]/"" '^-^ 

where the factor Xg^t comes from the path joining the two vertices labeled on the 
boundaries. Thanks to the Miermont bijection, this is also the generating function for 
our quadrangulations with a geodesic boundary of length 2{s + t) (up to a factor g^+^~'^ 
for the spurious added squares). Taking i = s + t, the formulas ( [A.l| ) and ( [A.2| ) thus 
enumerate the same objects, hence they must be equal. Equating (|A.1|) and ( |A.2|) with 
i = s + t yields 

^ _ [3]. [^ + 1L [^ + 1]. [^ + ^ + 3L ..ox 
[IL [« + 3L [t + 3l [s + t + 11 ^ ■ ' 



which is nothing but (1.5). 
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